Generalized minimum-uncertainty squeezed states 



E. ShchukinQ W. VogelQ and Th. Kiesefj] 

Arbeitsgruppe Quantenoptik, Institut fur Physik, Universitat Rostock, D-18051 Rostock, Germany 

Minimum-uncertainty squeezed states, related to a broad class of observables, are analyzed. Meth- 
ods for characterizing such states are developed, which are based on numerical solutions of ordinary 
differential equations. As typical examples we deal with nonlinear generalizations of quadrature 
squeezed states and deformed nonlinear squeezed states. In this manner one may derive those 
squeezed states which are directly related to given observables. This can be useful for optimized 
measurements at a reduced level of quantum-noise. 
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I. INTRODUCTION 

It is a well-known fact that the ground-state and the 
vacuum noise level of a harmonic oscillator and a mode 
of the radiation field, respectively, is required to fulfill 
the Heisenberg uncertainty relation. In fact, this noise 
level just defines the minimum uncertainty level of the 
quantum noise which is required to obey the uncertainty 
principle. It is important to note that the uncertainty 
principle sets a limit for the product of the variances of 
two observables, such as position and momentum of a 
harmonic oscillator and two orthogonal field quadratures 
for a radiation mode. Therefore it is a natural conclusion 
that one may find quantum states, for which the noise in 
one of the chosen pair of observables is reduced below 
the vacuum (or ground-state) noise level, at the expense 
of increased noise in the other observable [TJ [2] • Nowa- 
days such states are usually called squeezed states, in the 
early days of their study also the notion two-photon co- 
herent states was used [3[ [4] , since the structure of the 
unitary operator which leads to such state is formally a 
two-photon generalization of the coherent displacement 
operator. 

More than twenty years ago that the first successful ex- 
perimentell realization of squeezed radiation fields have 
been published [5] . There has been some interest in pos- 
sible applications of the noise reduction in a given ob- 
servable, for example in the context of intcrferometric 
detection of gravitational waves [6j [7] . It could be shown 
that the squeezing effect indeed improves interferometric 
[5J [5] and spectroscopic measurements 10J. Very recently 
squeezing could be realized with a reduction of the noise 
power by 10 dB [IT] , which makes the squeezed states 
indeed useful for gravitational wave astronomy. 

Based on this encouraging progress in the genera- 
tion and application of the so-called quadrature squeezed 
states, it is of some interest to rise the question of whether 
one may consider useful generalizations of the concept of 
squeezed states. For example, generalization of squccz- 
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ing has been proposed, which is based on the uncertainty 
relation of two general non-commuting observables |12j . 
The further investigation of squeezing in such a general 
sense is closely connected with the minimization of the 
uncertainty relation of two Hermitian operators. The 
problem of finding such generalized squeezed states re- 
quires to solve the minimization problem of the uncer- 
tainty relation for the two Hermitian operators defining 
the squeezing under consideration. The resulting mini- 
mum uncertainty states are then the generalized squeezed 
states. Generalized squeezing in this sense has been con- 
sidered for some special choices of the basic operators, 
such as amplitude-squared squeezing [13 and its higher- 
order generalizations [TH [15] , for a review see also [16] . 
Other generalizations of quadrature squeezing were based 
on the consideration of higher-order moments |17j . and 
on higher powers of the annihilation/creation operators 
in the squeeze operator [18]. 

In the present paper we will consider minimum- 
uncertainty squeezed states for two general noncommut- 
ing observables. Pure quantum states which fulfill this 
requirement will be constructed as the solutions of an 
eigenvalue problem, which can be analytically solved 
only in a few special cases, including quadrature squeez- 
ing and amplitude squared squeezing. For more general 
choices of the two Hermitian operators, we develop a 
systematic approach to find such states numerically in 
the Fock-Bargmann representation. This method is ap- 
plied to other types of squeezed states, such as general- 
ized quadrature squeezed states and deformed nonlinear 
squeezed states. It allows one to obtain and characterize 
the optimized squeezed states for a chosen observable. 
This may provide a powerful tool for optimizing a given 
measurement principle with respect to the relevant level 
of quantum noise. 

The paper is organized as follows. In Sec. 2 we provide 
some useful methodical details of the Fock-Bargman rep- 
resentation. The definition of the generalized squeezed 
states under study is introduced in Sec. 3. The gener- 
alization of quadrature squeezing is studied in Sec. 4, 
where in the quadrature operator the annihilation oper- 
ator is replaced by a function of the latter. In Sec. 5 
we reconsider the known effects of quadrature squeezing 
and amplitude-squared squeezing from the viewpoint of 
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our general method, in these cases we are able to pro- 
vide analytical solutions for the corresponding squeezed 
states. Deformed nonlinear squeezed states are intro- 
duced in Sec. 6, and some of their properties are ana- 
lyzed. A summary and some conclusions are given in 
Sec. 7. 



II. FOCK-BARGMANN REPRESENTATION 

In this section we discuss the Fock-Bargmann represen- 
tation of pure states and give some usefull expressions for 
the quantities we need below. Any pure quantum state 
= ^2n=o c n\ n ) can be written as the action of the op- 
erator ijj(aJ), which is a function of the creation operator 
only, on the vacuum state 



1^=^)10), 
where the function ip(z) is defined via 

+oo _ _ 
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This series converges producing an entire analytical func- 
tion. The representation of quantum states by means of 
entire analytical functions according to Eq. (fTl) is referred 
to as Fock-Bargmann representation [T^J I20] 7 The scalar 
product {ipi\ip2) m this representation reads as 



(^ 2 >= / ri(z)Mz)e- lzl2 d 2 z 
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in particular, the normalization of a state is given by 
the following condition: 



(V|V)= / |^)| 2 e- |z|2 rf 2 z=£- 

J n=0 



|^ (n) (0)| s 



= 1. (4) 



If an entire analytical function ip(z) satisfies a weaker 
condition 



\^(z)\ 2 e-^ 2 d 2 z < + 



~Xj. 



(5) 



then the state defined by Eq. ([I]) is not normalized, and 
its normalization is the number M defined by Eq. (|5| . 

A function g(a') of the creation operator is the mul- 
tiplication by g(z) in the Fock-Bargmann representation. 
Due to the equality 



(6) 



the annihilation operator a corresponds to the derivative 
d/dz. The relation (All from |A| shows that 
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e a ^)\0) = £ ¥ ^ W (a+)|0) = ^ + a)|0), (7) 
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thereby the operator e aa is the shift of the argument in 
the Fock-Bargmann representation. The relation 



M >(a t )|0)=^ T )|0) 



(8) 



shows that the operator /1" corresponds to the scaling 
of the argument. The moments (d^ n a m ) in the Fock- 
Bargmann representation are readily calculated as fol- 
lows: 



k=Q 



(9) 



Usually it is much easier to use the discrete versions of 
Eqs. ([3]), Q and ^ than to calculate the corresponding 
double integrals. Using Eq. ([7]) one can get the following 
expression for the normally-ordered characteristic func- 
tion $(/?,/?*) = (e^ ' e-P* h ) of the state (0: 



$(/?,/3*) = / ijj*( z + f3)i)(z- /3*)e"l z l 2 d 2 z 



g T/^ fe )(/3)^ (fc) (-/3*) 

k=0 



(10) 
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Note that this expression is in agreement with Eq. ^ 
since the moments (a^ n d m ) can be expressed in terms of 
the characteristic function <&(/?,/?*) as follows: 

<9™+ m $(/3,/?*) 



(a tn a m ) 



d/3 n d(-/3*y 



(11) 
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Below we need to calculate the moments of the squeezed 
variant of the state ([I]) 



\tp sq ) = s(om = s(OH^)\o). 



(12) 



where S(£) 



(Ca 2 -ca+ 2 )/2 



is the squeezing operator. 



The moments of this state can be obtained from Eq. (Ill 



and the following general relation for the characteristic 
function <f> sq (/3, j3*) of the squeezed variant of any quan- 
tum state: 



$ sq (/3,/3*) = $(/?', /3 '*) e -kl 2 l/5| 2 -pRc(^/3 2 ) j 



(13) 



where j3' = (i[3+isf3* . In some examples considered below 
there is no closed analytical expression for the character- 
istic function (to our knowledge), but in all examples it 
is possible to get such expressions for the moments using 
Eq. (JoJ) (though sometimes these expressions are really 
huge). 



III. GENERALIZED SQUEEZING 

Let us now consider the definition of generalized 
squeezed states in more detail. The calculation of the 
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properties of such states will be reduced to the problem 
of solving ordinary differential equations. Some examples 
will be studied in the next section. 

The general uncertainty relation for two Hermitian op- 
erators F and G of equal dimension reads as follows: 



AFAG>±\([F,G})\, 



(14) 



where AA = {(A A) 2 ) 1 / 2 is the dispersion of A. In this 
work we study the states which minimize this uncertainty 
relation, i.e. the states which satisfy the equality 



AFAG=-\{[F,G]}\. 



(15) 



Unless AF — AG, exactly one of the following inequali- 
ties is valid: 



According to Eq. ([I]) we can look for the solution (a\ip) 
in the form 



*\ p -\a\ 2 /2 



(21) 



whereby Eq. ( 20 I can be simplified to an ordinary differ- 



ential equation of a complex variable 



F 



d 

da* 



+ iXG a 



d 

da* 



jjj(a*) = /3i)(a*). 



To get an ordinary differential equation of a real variable, 
let us represent a* in polar coordinates as a* — re~ %v . 
The unknown function ip( a *) can t* e considered as a func- 
tion of the radius r for a fixed phase <p> 



</>(a*) = ^(re-*) = tf> v (r). 



(23) 



(AF) 2 <i|([F,G])| or (AG) 2 < ±\([F, G])\. (16) 

A state that satisfies any of these inequalities was called 
generalized squeezed state f [T21 12"T] ). 

In this work we deal with pure states only. Any solu- 
tion |?/>) of the eigenvalue problem 



(F + iXG)\ip) = (3\ip), 



(17) 



The derivative V^( r ) (with respect to r) can be calculated 
with the help of the standard chain rule 



dip(a*) da* _ d^(a*) _ 4 



da* dr 



da* 



(24) 



thereby the derivative with respect to complex argument 
a* is related to the derivative with respect to the radius 
(for the phase fixed) via 



where A is a positive real number and (3 is arbitrary com- 
plex, also satisfies Eq. (15 1 cf. (32]. The relation between 



the dispersions AF and AG for such a state reads as 



AF = AAG, 



(18) 



da* dr 



(25) 



Now Eq. ( 20 ) can be written as an ordinary differential 
equation 



so the parameter A plays the role of the degree of squeez- 
ing. For < A < 1 the first of the inequalities ( 16 I is 
satisfied, and for A > 1 the second one is. The solutions 



of Eq. (17 1 for A = 1 are unsqueezed (in the generalized 
sense under consideration). Nevertheless the resulting 
state can be a nonclassical one. 



In general, it is impossible to solve Eq. (17) analyti 



cally, but one can rewrite it as an ordinary differential 
equation, whereby making it possible to solve it numer- 



ically. To transform Eq. (17 1 to an ordinary differential 



equation note that for a coherent state \a) we have the 
relations 



(a\a 



a (a\ 



{a\a 



d 

da* 



(a 



(19) 



We assume that the operators F and G are written in 



the normally-ordered form. Using the relations (19 1 we 



can get the following differential equation for the scalar 
product {a\ip): 



_d_ 

da* 



+ i\G 



' 2 



d 

da* 



(20) 



F ( re- iv ,e iv ^\ + iXG ( re - ** e iv ~ 
\ dr J \ dr 

= /?Vv( r ), 



ipjr) 



(26) 



or, more precisely, as a family of equations parameter- 
ized with the phase tp G [0,27r]. Solving this equation for 
all phases and combining the solutions ip<pi r ) together we 
get a solution tp( a *) °f Eq. (22 1. The initial conditions 



to Eq. (26) depend on the phase tp, but the final solu- 
tion ip(a*) must be an analytical function of a* , thereby 
one cannot take arbitrary functions of ip as initial condi- 
tions to Eq. (26 1. For ^(a*) to be analytical the initial 
conditions must be chosen as follows: 



Vi fc) (o) = Cfce -^, 



fc = 0,l,.. 



(27) 



with Cfc being arbitrary complex numbers. The deficiency 
of this method is that, in general, for arbitrarily chosen 
initial conditions (27 1, the solution ip(a*) we get in this 



way is not normalized, one must numerically take the in- 
tegral ([5]) to find the normalization constant. In the fig- 
ures of the numerically calculated Q-functions presented 
below, we did not normalize the corresponding states, 
since in those figures it is just a matter of scaling, which 
is unimportant for our purposes. 
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IV. NONLINEAR QUADRATURE SQUEEZING and the corresponding differential equation (26 1 as 



In this section we introduce a nonlinear generalization 
of quadrature squeezed states and we prove their nonclas- 
sicality. First we define the nonlinear generalization of 
the quadratures to be studied. The known special cases 
of quadrature squeezing and amplitude-squared squeez- 
ing will be reconsidered from this generalized point of 
view in the following section. 

Any pair of Hermitian operators F and G can be rep- 
resented in the form 



F=f + P, G=-i(f-p) 



(28) 



for some operator /, one may just set / = (F + iG)/2. In 
all the examples considered below we use this representa- 
tion with different choices of /. Note that the operators 
that are defined in the following may be considered as 
a direct nonlinear generalization of the quadrature op- 
erators, which are recovered in the linear special case, 
f = a. 

Let us consider now the case of the operator / being 
a function of the annihilation operator only, / = /(a). 
In order to not overload the notation and without loss 
of generality, we assume in the following the function 
f(z) to be real (i.e. for real argument z its value f(z) is 
also real), which implies that the relation f(a)^ = /(a^) 
is valid. Let us write the generalized squeezing condi- 
tions given by the inequalities (16) explicitly The first 
inequality is 



(AF) 2 <i|([F,G])|. 



(29) 



The variance (AF) 2 reads as 



(AF) = ((AfY + (A/) + A/' A/ + A/A/'), (30) 

and for the commutator of F and G we have the equality 
([F, G]) = 2i([f, /t]). In|X]it is shown that that the lat- 
ter commutator is always nonnegative (in the case under 
study): ([/,/t]) = ([A/, A/t]) > 0. Now the inequality 
(29 1 can be written as follows: 



((A/t) 2 ) + ((A/) 2 ) + 2(A/tA/) < 0. (31) 

Since / — /(a), the left-hand side of this inequality is 
just (: (AF) 2 :), so the condition (29) reads as 



(AF) 2 :) < 0. 



(32) 



From the second inequality of (16 1 we can conclude in 
the same way that (: (AG) 2 :) < 0. We see that gener- 
alized quadrature squeezing always implies nonclassical 
behavior. 

The eigenvalue problem (JTtJ) reads as 

((1 + X)f(a) + (1 - A)/(ot)) |V) = W), (33) 



dr ) ^ 



A - 1 
A + l 



fire-**) + 







A+l 



(34) 

In general, for an arbitrary function f(z), it is impossible 
to solve these equations analytically. But in a special 
case of A = 1 it is possible to find the general solution, 
provided that we know how to find the roots of entire 
functions. In the case of A = 1 Eq. (33 1 simply reads 
as f(a)\ip) = 7IV'), where 7 = (3/2. Each root a to the 
equation f(a) — 7 gives a partial solution of the form 
Ffc a _i(at — a) I a), where k a is the multiplicity of the root 
a and Pk a -\{z) is an arbitrary polynomial of the degree 
k a — 1. In particular, any simple root (of multiplicity 1) 
gives the solution which is proportional to the coherent 
state I a). The general solution of Eq. (33 1 for A = 1 
is a linear combination of all the partial solutions 



/(")=7 



L (a f - a*)\a), 



(35) 



where the sum here is taken over all the roots of the equa- 
tion f(a) = 7. If all the roots of this equation are simple 
then the general solution ( |35| is just a linear combina- 
tion of coherent states. Below we consider only polyno- 
mial functions f(z). Note that, though there is no gen- 
eral analytical expression for the roots of a polynomial 
of a degree greater than four, it is possible to find out 
whether all the roots of the polynomial in question are 
simple or not (a polynomial whose all roots are simple 
is called separable), which answers the question whether 
the sum ( [35] ) contains only coherent states or not. A 
polynomial is separable if and only if a special determi- 
nant constructed from the coefficients of the polynomial 
is not equal to zero. The details see, for example, in [25] . 
It is clear to see that Eq. ( 34 ) with a first-order order 



polynomial /(z) corresponds to the definition of quadra- 
ture squeezing, and it is easy to prove that a second-order 
polynomial leads to the definition of amplitude-squared 
squeezing. That is why it makes sense to consider poly- 
nomials of third and higher orders. As an example, let us 
consider the polynomial f(z) = z 3 + z. The differential 
equation ( 34 ) with this polynomial reads as 



A+l 



(r 3 



3i(p 



re 



-up 



) + 







A + l 



(36) 



The contour plot of the Q-function of a solution of this 
equation is shown in Fig. [I] for some values of the squeez- 
ing parameter A. 

In the next sections we consider the two cases (quadra- 
ture squeezing and amplitude-squared squeezing) in a 
more general context, allowing the squeezing parameter 
A to be complex. We find the condition which guarantees 
that solutions of the corresponding equations are normal- 
izable (and thus can be identified with physical quantum 
states) and show that only for real A are these states 
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(a)A = 0.1 (b)A = 0.2 (c)A = 0.5 



(d)A = 1 



(e)A = 2 



(f)A = 5 



(g)A = 10 



FIG. 1: The Q-function corresponding to the solution of Eq. ( 36 1 with the initial conditions V (0) = 1, 9' (0) = 0, 6Z(0) = e 



minimum uncertainty. The details of the calculation are 
given in the appendices. 



V. QUADRATURE AND 
AMPLITUDE-SQUARED SQUEEZING 
REVISITED 



The condition of convergence \z\ < 1 is equivalent to the 



positivity of ReA, thereby the state (38 1 is normalizablc 



if and only if ReA > 0. The normalization M is given by 



TV 2 



Re 



|A + 1| 



exp 



A*- l 
A+ 



4ReA 



(41) 



Let us now reconsider the well-known cases of quadra- 
ture squeezing and amplitude-squared squeezing in the 
framework of our approach. For these cases we can solve 
the problems analytically. In the next section we will 
further develop our approach for the study of deformed 
nonlinear squeezing. 



A. Quadrature squeezing 

In the case of quadrature squeezing the operator / is 
simply the annihilation operator f = a. The operators 
F and G defined by Eq. (28 1 are then two orthogonal 



quadratures: F — x and G = p. The differential equation 
(22 1 in this case reads as 



da* 



A-l 
ATT* 



(3 



A + l 



(37) 



This equation can be easily solved which gives us the 
eigenstates of the corresponding eigenvalue problem ( 33 ) 



= N exp 



1 A- 1. 

2 A+ l' 



t2 



A + l 



(38) 



To calculate the normalization Af using the relation ( Al ) 
lg equa 

(z/2)» 



we need the following equality for Hcrmitc polynomials: 



1 



n=0 



VT^ 2 



( 2xyz — (x 2 + y 2 )z 2 



(39) 



The series on the left hand side converges only if \z\ < 1. 
Here for x, y and z we get 



/9 



x 



A+l 



(40) 



For A = 1 the state fl38| is just a coherent state |/3/2). 

For real A the state ( |38| minimizes the uncertainty re- 
lation ( 14 1: AxAp = 1. Let us see what happens when 



A is complex. To find the dispersions of x and p for the 
state ( 38 ) we need to calculate its normally ordered mo- 



ments up to second order. Note that it is straightforward 
to calculate the antinormally ordered moments {a n a) m \ 



for the state (38 1 due to the following relation: 



(o n o ,m ) =A/" 2 (A* + 1)"(A + 1)' 



d(3* n d(5 m 



(42) 



It is possible to express these moments in terms of 
Hcrmitc polynomials and then obtain normally ordered 
moments in general, but we need only second-order 
moments, which can be easily obtained directly from 
Eq. ( 42 ) . For the simplest moment (a) we have 



Re(A/3*) + flm/3 
2ReA ' 



(43) 



The expressions for (a ) and (cva) are more lengthy, so 
we present here only the final expressions for the disper- 
sions 



(Ax) 2 



ReA' 



(Ap) 2 = 



1 

ReA' 



(44) 



To see how much the product of the dispersions differs 
from the lowest possible value 1 let us calculate the dif- 
ference 



(Ax) 2 (A P ) 2 



|A| 2 
Re 2 A 



1 = tan 2 



(45) 



where <& = argA. We see that the state (38) is minimum 



uncertainty if and only if A is real, and the larger $ is, 
the stronger the equality AxAp = 1 is violated. This 
is closely related to the problem of optimal choice of the 
phase of the quadratures adjusted to the principal axes of 
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the squeezing ellipses, cf . [I] . When A tends to the imag- 
inary axis (i.e. when ip tends to ±7r/2), the difference 
( 45 I tends to infinity. 



It is also interesting to see how the squeezing of the 
state (38 1 depends on A. To see this we must calculate 
the minimum min ¥ ,(: (Ax v ) 2 :) of the general quadrature 
-i<p + at e **>. it i s not difficult to do and after all 



ae 



simplifications this minimum reads as follows: 



min(: (Ax v ) 2 



2IA-1I 



IA + 1I + IA-1I 



(46) 



We see that the state under study is always squeezed, 
except the case of A = 1, when it is just a coherent state, 
as it has already been mentioned above. 



B. Amplitude-squared squeezing 

In the case of amplitude squared squeezing the opera- 
tor / reads as / = a 2 . The differential equation (22 1 is 
now a second-order equation 



d 2 ip(a* 
da* 



A - 1 
A + l 



cr 2 + -r^— I ip(a*) = 0. (47) 



A + l 



It has the following two linearly independent solutions 
(called even and odd, for self-evident reason): 



if) e (b; c; a* 



> \F X ( 6; \;2ca* 2 



ipo(b; c; a*) = a* e 



1 3 

i-Fi I b + -; -;2ca 
' 2 2 



(48) 



where iFi(b; c; z) is the Kummcr function denned via 



1 F 1 (b;c;z) = J2 



k=0 



(C)fc k\ ' 



(49) 



with (6)^ being the Pochhammer symbol expressed in 
terms of gamma functions as (b)k = T(b + k)/T(b) = 
6(6 + 1) . . . (6 + k — 1). The parameters b and c read as 



A- 1 
ATI' 



fe =4 (1 







vx 2 



(50) 



Note that there are two choices for the parameter c which 
differ in sign only (since each nonzero complex number 
has two square roots). When the sign of the square 
root -\/(A — 1)/(A + 1) is chosen, the square root y/X 2 — 1 
for the parameter b must be calculated according to the 
equality \/A 2 - 1 = (A - 1)/a/(A - 1)/(A + 1). liMit is 
shown that both the solutions ip c and tjj , Eq. ( |48[ ), do 
not depend on the choice of the sign of the square root 
for c. 

It is straightforward to write the quantum states which 
correspond to the functions (48) in the Fock-Bargmann 



representation, but for us it is more convenient to repre- 
sent them in another form. In fact, as shown in[Bl the 



even and odd solutions of the eigenvalue problem ( 33 I 
read as 



\^ c ) = KS(0iFi (b;\;va^ |0), 
\i>o)=KS(0 1 F 1 U+I;| ;w a+ 2 ) |1>, 
where the parameters £ and v are defined via 



(51) 



e^tanhiei^/^j, v = 



1 



A-l 



A+l 



(52) 



and the normalizations J\f c and J\f Q are as follows: 

1 



AC 2 
AC 2 



2-Fi I b,b';-; 



4« 



i/ 7 i ! b+\,b* 



1 3 

2' 2' 



■Av 



(53) 



Note that for real A the states (|51| are exactly the ones 
obtained in [24]. The contour plot of the Q- function of 
the even state is shown in Fig. [2j 

In the case of quadrature squeezing there arc relatively 
compact analytical expressions for the general normally 
ordered moments of the solutions of the corresponding 
eigenvalue problem ( 33 ) . In the present case the situa- 



tion is much more complicated. To our knowledge, there 
is no analytical expression for the characteristic function 
of the states (51), but using Eqs. (Ill and ^ together 
with Eq. (13) it is possible to find analytical expressions 



for normally-ordered moments of these states. But these 
expressions are really huge, for example, only the single 
moment (a 4 ) (which is needed to calculate the disper- 
sions of F and G) expands to a whopping whole-page 
expression. Let us present here only the mean photon 
number of the even state 



= 4|6| 



A 2 II f\\+l\-\\-l\ 



ReA 
2 Fi (6+1,6* 



lA+ll + IA-ll 



l;f;4|c|- 



2 F 1 (6,6*;i;4|c| 2 ) 



(54) 



JA 2 



II 



ReA 



^Re6- 



|A-1| 



lA+ll-lA-ll 



Note again, that for real A this coincides with the expres- 
sion for the mean photon number obtained in |24j . We 
found that the states (51 ) are minimum uncertainty only 



for real A, and the larger the phase of A the stronger is 
the violation of the corresponding equation (15), though 



in this case the strength of the violation does not depend 
only on the phase of A, as it was in the case of quadra- 
ture squeezing. Calculating the moments with Mathe- 
matica we experienced dramatic loss of accuracy work- 
ing with machine precision, so it was necessary to use 
high- precision numbers to get the correct results. 
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(a)A = 0.1+i (b)A = 0.2 + j (c)A = 0.5 + i (d)A = l + j (e)A = 2 + i (f)A = 5 + i (g)A = 10 + i 



FIG. 2: The Q-function of the even state (51 1 for (3 = 5 and different A 



VI. DEFORMED NONLINEAR SQUEEZED 
STATES 

Here we consider a more complex case of the operator 
/ from Eq. (28 1 to be of the form / = g(fi)a, where we 



again assume g(z) to be real. The eigenvalue problem 
( 17 1 in this case reads as 



(1 + X)g(n)a + (1 - \)tfg(h)) = /3|V>- (55) 



To formulate the corresponding differential equation ( 26 ) 



we need to normally order the functions of the photon 
number operator. In Owe prove the following relation: 



9{n) = £ 



fc=0 



(A fc g)(0) 
fc! 



(56) 



where A is the difference operator defined via (Ag)(n) = 
g(n + 1) — g(n). The powers of this operator calculated 
at zero read as 



(A fe 5 )(o) = ]T(-i; 



k—i 



i=0 



<?(*)• 



(57) 



Now we can write the differential equation ( 26 1 explicitly 

+00 



fe=0 



d k ip v (r) 
dr k 



(58) 



where the functions C^(r) = Cfc(A,r, ip), k > 1, are de- 
fined via 



C„(r) = (1 - X)g(0)re-^, 

(A k ~ 1 g)(0) t_ x 



C fc (r) = (1 + A) 



+ (1-A) 



(fc-1)! 
(A fe 5 )(0) 



(59) 



fc! 



The initial conditions are given as in Eq. (27 1, but in 



this case the coefficients of the equation are not constants 
and that is why the initial conditions cannot be chosen 
arbitrarily. It is shown in [D] that there arc the following 
relations for the initial conditions to Eq. (|58|), k > 0: 



(60) 



fc(l - X)g(k ~ l)e-^^ fc - 1 )(0) - Wf(0) 
+ (1 + A).g(fc)e^ +1) (0)=0. 



Usually (but not always, as we will see shortly) this 
means that the solution of Eq. (58 1 is unique. 

In the very special case of A = 1 the solution of the 
equation (55) reads as (see |2"5"] ) 




7" 



\ Vn\ 5(0) • • • g(n 



1) 



(61) 



provided that all the numbers g(n), n = 0, 1, . . . are not 
equal to zero. As an example let us take g(h) = n. In 
this case Eq. ( 55 1 reads as 



(1 + A) re' 



dr 2 



+ (l-A)r 2 e -^^^=/?^(r). 

The constraints ( 60 1 on the initial conditions in this 
case simply read as /?Vv(0) = 0j thereby if (3 ^ then 
Vv(0) = and if [3 = then both the derivatives W(0) 
and ?/^(0) can be arbitrary. Since it is a second-order 
equation and if one of its initial conditions is fixed, then 
its (normalized) solution is unique. The Q- function of 
the solution is shown in Fig. [3j In the case of A = 1 it 
is possible to solve Eq. d62| analytically, but we cannot 



directly use the expression given in Eq. (61 1 since in this 
case g(0) = 0. If (3 ^ then the solution for A = 1 reads 
as follows: 

= oF 2 (;l,2; |/?| 2 /4r 1/2 Fi(;2;/3aV2)|l), (63) 

where oFi(;b;z) is the regularized confluent hypergeo- 
metric function defined via 

g, , , » oFi(;b;z) ^ I z k 

F 1 (;b;z) = — -j-^ — , Fx(;b; z) = ^ 



T{b) 



k=0 



If (3 — then the solution is just a linear combination of 
the first two Fock states \ip) = cq\0) + Ci|l), hence the 
solution is not unique. 



VII. SUMMARY AND CONCLUSIONS 

In the present paper we have studied generalized 
squeezed states, which are based on the minimum of the 
uncertainty relation for two general Hermitian operators. 
The problem of finding the related generalized squeezed 
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FIG. 3: The Q-mnction corresponding to the solution of Eq. (621 with initial conditions 9 V (Q) = 0, 0' V (Q) = e 



states has been reduced to solving ordinary differential 
equations which are obtained in the Fock-Bargmann rep- 
resentation. These equations can be solved analytically 
only in some special cases, such as quadrature squeez- 
ing and amplitude-squared squeezing. For more general 
cases we have developed techniques to solve the differen- 
tial equations numerically. 

To illustrate the power of our approach, we have stud- 
ied two examples of generalized squeezed states of types, 
which to our best knowledge have not been considered 
so far. The first type is a nonlinear generalization of 
the quadrature squeezed states, where in the definition 
of the quadrature operators the annihilation operator is 
replaced by a function of the latter. The second type is 
a deformed nonlinear squeezed state, which is defined on 
the basis of the quadratures of a deformed algebra. For 
both types of generalized squeezed states we illustrate 
their properties by calculating the phase-space distribu- 
tions, i.e. the Q-functions. 

In conclusion, we have studied several types of general- 
ized squeezed states. The reduction of the quantum noise 
level of such states is related to different types of Hermi- 
tian operators. In this sense such kinds of squeezed states 
may become interesting when one wants to improve spe- 
cial measurement schemes. In such cases it may be of 
interest to find the best squeezed states in relation to 
the observable to be detected by a given device. When 
knowing the measurement scheme, our method is useful 
to characterize the squeezed states which are adjusted to 
the observation scheme. This may be a first step towards 
the preparation of these states and their applications for 
optimized measurements at a reduced level of quantum- 
noise. 



APPENDIX A: NORMAL ORDERING 



In this appendix we prove the following relation: 



to get the equality 

n 

a"/(6 t )=E(l t )/ (fc) (« t )a n - fc , (A2) 



k=0 



by ind uctio n. Then we can calculate the left-hand side 
of Eq. ( Al |, assuming g(x) = J2n=o 9nX n 



+oo n 



n=Ok=0 ^ ' 

+oo 1 +oo . 



n—k 



(A3) 



n—k 



= E^/ (fc W fe) (a), 



k=0 



which gives us the desired result. Here we used the fol- 
lowing symbolical rule for double sums: 



+oo n +00 +oo 

EE-EE- 

n=0 fc=0 k=On=k 



(A4) 



If the function /(z) is real then the commutator 
[f(a),f(a')} reads as 

[/(a),/(at)] = ^-/( fe )(a)t / W(a), (A5) 



fc=i 



so that it is nonnegative. 



APPENDIX B: AMPLITUDE SQUARED 
SQUEEZING 

In this appendix we prove that the states given in the 
Fock-Bargmann representation ( 48 1 can be rewritten in 



the form of Eq. (51 1. First of all, we show that the func- 



tions in Eq. ( 48 1 do not depend on the choice of the sign 



g(a)f(a ! ) = ~p k '(a^)g^(&), (Al) for the parameter c in Eq. (501. There are two possible 



fc=0 



where f(z) and g(z) are entire functions. From the ubiq- 
uitous bosonic commutation relation [a, al] — 1 it is easy 



values for c, which differ in sign, denoted here as c+ and 
c_ so that c + = — c_. To c + and c_ correspond two 
values for b, b + and 6_ respectively. As one can easily 



see from Eq. (50) they are related as b + = 1/2 — 6_. It 
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is well known that the Kummer function 1F1 satisfies to 
the following relation: 

1 F 1 (b;c;-z) = e~ z 1 F 1 (c-b;c;z). (Bl) 

Using this relation we can write 

V> c (b-; c_; z) = e c + z \Fi (b-, i; -2c+z : 

=e -^ lFl (i. b _.,i, 2c+ ^ m 

= ^c(b+;c+;z). 

The same is also true for the odd function ip (b; c; z). 
Now we prove the following relation: 



x e 



va^ 2 



2 2 ^* w ) 
where the parameters /i, and £ read as 



(B3) 



|0>, 



At = cosh|f|, ^ = e iars «sinh|£|, C = -• (B4) 

A* 

Note that the squeezing operator can be written as 
S(0 = e? where the operators K± and Kq are 

defined to be 



5t2 



a 2 



These operators satisfy to the following commutation re- 
lations: 

[k , k±] = ±k ± , [k-,k+] = 2k , (B6) 

and these three operators generate the Lie algebra of the 
group SU(1, 1). The normal form of the squeezing oper- 
ator is given by (cf. [26] ) 



S(0 = e 



<K+ 



2K 



oCK- 



1 / Ca f2 C*a 2 M-l. 
— : exp 1 1 n 



(B7) 



We will transform the left-hand side of (B3) step by 



step. First, we must calculate the expression \ip) — 
e C"a / 2 1 _F 1 (6; 1/2; va^ 2 ) |0) and this is the only nontrivial 
step in the whole process. According to the relation ( Al ) 
we can write it as 



n=0 



dz r 



The derivative of the exponent at zero reads as follows: 



d n e C Z 2 /2 



dz n 



z=0 





(C*/2) 



n/2_ 



WW 



if n is odd 
if n is even, 



(B9) 



thereby only terms with even n survive in the sum ( B8 1 



An even order derivative of the Kummer function reads 
as 



-■vz z \ = (4<(6) fe 1 F 1 ib + k;-;vz- 



(B10) 

Now we can further proceed with the expression for 



+ 00 

£(% 

k=0 



(2(*v) k 



1 



ita 



(Bll) 



According to the definition ( 49 ) of the Kummer function 
we have the equality 



(B12) 



and upon substituting it into the previous relation we 
finally get 

{ ~ )= g (2Cv) k (b) k+m (v^ 2 r^ 



k,m—0 



k\ (l/2) m ml 



1 +°° / 

(i -2c*w) & ^ It 

v s I m=0 



*t2 



(b)r 



2(*vJ (l/2) r 



1 f, 1 va 

Fx | & 



(1 -2C*v) b 



t2 



2 1- 2C*v 



|0). 



(B13) 



Here we used the relation (&)&(&+ fc) m = (6)fe+ m and the 
following equality: 



£(&)*H 



fc=0 



fc! 



(1-x) 



b-\-m ' 



(B14) 



We have just obtained the following expression for 



I* 



1 ,« U 1 



(1 -2Cv) b 



2' 1 - 2C,*v 



(B15) 



Then we must apply the operator (l//i)™ +1 / 2 to \ip). 
As we already mentioned, this is equivalent to scaling 
the argument of the corresponding function in the Fock- 
Bargmann representation. Applying this scaling to the 
left-hand side of the previous equati on a nd multiplying 
by e - ^ a+2 / 2 , we finally arrive to Eq. (B3). 



Now we must find £ and v such that the right-hand side 



of Eq. (B3 1 is exactly the even state defined in Eq. (51 1 



|0) (B8) K is easy to see that £ and v defined in Eq. ( 52 1 have this 



property. This finishes the proof of the representation of 
even and odd states in the Fock-Bargmann representa- 
tion given by Eq. (48 1 in the form of Eq. (51 1. 



Now we calculate the normalization of the even and 
odd states under study. It is much easier to do when 



these states are represented in the form of Eq. (51 1, for 



10 



this reason we have transformed the states into this form. 
The normalization J\f c of the even state can be calculated 
as 



+00 1 , . 



n=0 



(n) 



y . , x Fx ( b; --vz 2 

I z=0 \ z / z=0 



1 



(n) 



The derivatives in this expression read as 



(B16) 



if n is odd 



1 2 \ (n) [ 
1/1 {b; 2' Vz2 ) z=0 ^\(4vr^(b) n/2 ifniseven, 



(B17) 

thereby the normalization Af c can be written as follows: 



+00 



« (i6i«i 2 r 



(2m)! 



(B18) 



Using the relation (2m)! = 4 m m! (1/2) , we finally get 
Eq. ( 53 1 . In the same way one can obtain the normaliza- 
tion A/q of the odd state. 



APPENDIX C: NORMAL ORDERING OF 
PHOTON-NUMBER OPERATOR FUNCTIONS 

In this appendix we prove the following equality: 

g{h) =: (e AA g)(0) :, (CI) 

where f(z) is an entire function. Let us start with the ex- 
pression for normally-ordered powers of the photon num- 
ber operator 

m 

: n m := n(n - 1) . . . (h - m + 1) = £ s(m, k)h k , (C2) 

k=0 

where s(m, k) are the signed Stirling numbers of the first 
kind. One can invert this relation and express ordinary 
powers of h in terms of normally-ordered ones 



(C3) 



fe=0 



where S(m, k) are the Stirling numbers of the second 
kind. Explicitly these numbers read as 



(C4) 



This is the number of ways to partition a set of m el- 
ements into k non-empty subsets and hence it is equal 
to zero if m < k (see [27]). Using the Taylor expansion 
g(z) = X)m=o 9mZ m of the function f(z) and replacing z 
with h we can write the operator g(h) as 

+00 +00 m 

9{n) = £ g m n m = £ 5m £ S ( m > k ) '■ ^ '■ 

m=0 m=0 k=0 

+OO / +00 \ +00 

= E E )■''■ E F * : : • 



fc— \rn—k 



k=0 



(C5) 



Now we have to calculate i 7 ^ defined via 

+00 k s,\ +00 

F k =J2 S{m, k)g m = - ]T(-l) fe -* . £ 9mf 



-k i—0 m—k 



(C6) 



The second sum can be represented as follows: 

+00 k—l 

E 9^ m = g(i) - £ g m i m , (C7) 



m=0 



so that .Fi. can now be written as 



Fk ^E^^ftM ~ E (C8) 

i=0 ^ ' m=0 

It has been mentioned above that S(m, k) = for m < k, 
so the last sum in this expression is zero, and finally we 
have 

F » = hl2(- 1 ) k ~ i ( 1 f)9(i} = ^(A fe 9 )(0), (C9) 



which completes the proof of Eq. (CI I 



APPENDIX D: THE RELATIONS FOR THE 
INITIAL CONDITIONS 



In this appendix we prove the relation ( 60 I for the 



initial conditions for Eq. ( 58 1 . It is easy to see that 



the derivatives at zero of the functions Cfc(r) defined in 
Eq. (59) read as 



(l + A)(A fc - 1 . 9 )(0)e^ 



= k- 1 



Cr (°)= { (fc + l)(l-A)(A fe «?)(0)e-^ m=k+\ 
m^k±l 

(Dl) 



Differentiating Eq. (58) m times at r = we get the 
following relation: 



+00 m 



k=0 J=0 



EE ( 7 ) ci j) m¥ +m - j Ho) = /?^ m) (o). (D2) 



Using the expressions (Dl ) this relation can be simplified 
as 

rn , . 

E ? (^ ) 1 (oX +1 )(o) + cji ) 1 (o)^ m - 1 )(o)) 
= M m) (o). 

(D3) 

The coefficient in front of the derivative ip^ 1 (0) can 
be further simplified as follows: 

m / \ m f \ 

E ? jCi+iW = (!+ A )^E ? ( AJ 5)(o) 

3=0 ^ ' i=o ^ 7 ' 

= (1 + A) e ^((l + A) m 5 )(0) 

- (1 + X)e^(E m g)(0) = (1 + A) 3 (m)e^, 

(D4) 
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where E — 1 + A is the step operator which acts as 
(Eg)(n) = g(n + 1). In the same way one can calculate 

the coefficient in front of i[)^, n ^(0) and get 



This completes the proof of the relation ( 60 ) 



E 

j=0 



CR(0) = m(l - X)g(m - l)e 



•lip 



(D5) 
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